arXiv:1503.01932vl [cond-mat.stat-mech] 6 Mar 2015 


Incommensurate Single-Angle Spiral Orderings of Classical Heisenberg Spins on 

Zigzag Ladder Lattices 

Yu. I. Dublenych 

Institute for Condensed Matter Physics, National Academy of Sciences of Ukraine, 1 Svientsitskii Street, 7 9011 Lviv, Ukraine 

(Dated: March 9, 2015) 

Exact and rigorous solutions of the ground-state problem for the classical Heisenberg model with 
nearest-neighbor interactions on two- and three-dimensional lattices composed of zigzag (triangular) 
ladders are obtained in a very simple way, with the use of a cluster method. It is shown how the 
geometrical frustration due to the presence of triangles as structural units leads to the emergence 
of incommensurate spiral orderings and their collinear limits. Interestingly, these orderings are 
determined by a single angle (along with the signs of the interactions between neighboring spins); 
therefore, the term “single-angle spiral orderings” is proposed. 
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In geometrically frustrated systems, the energies of all 
interactions cannot be simultaneously minimized because 
of geometric constraints [lj. There is no general algo¬ 
rithm for determining ground states even for classical 
frustrated systems. For this reason, the ground-state 
problem for many classical Heisenberg spin models still 
remains a challenging task for theorists [§, 3|. The most 
known and used method for these purposes is that of Lut- 
tinger and Tisza [3] and its generalizations [lj. However, 
it appears not to be useful for complicated Hamiltoni¬ 
ans (for instance, with an external field or a biquadratic 
exchange) and difficult to apply to non-Bravais lattices. 
There is also a cluster method proposed by Lyons and 
Kaplan half a century ago Q, but it is “rather unknown,” 
although simple and intuitively clear [3, :8). We developed 
this method in our studies of the ground states for some 
Ising-type models [MU- Here, we consider the classical 
Heisenberg model on a set of lattices for which, despite 
the presence of frustration and a non-Bravais character 
of the lattices, a solution of the ground-state problem 
can be obtained in a remarkably simple and clear way, 
with the use of the cluster method. These are two- and 
three-dimensional lattices composed of zigzag (triangu¬ 
lar) ladders. There are quite a few compounds where 
magnetic atoms are arranged in lattices of this type 0 - 
IH. Making use of the cluster method, we clearly and 
rigorously show how the presence of structural triangles 
and, hence, the geometric frustration, leads to the emer¬ 
gence of incommensurate spiral orderings. An important 
point is that, despite the presence of several interaction 
parameters, these orderings are determined by a single 
angle (together with the signs of the interactions between 
neighboring spins along the ladder rungs). Therefore, we 
introduce the term “single-angle spiral orderings.” 

Since triangular ladders are composed of identical 
isosceles triangles, let us start with a single isosceles tri¬ 
angle and classical spins (unit 3 vectors) at its vertexes 
(Fig. 1). L and K are the strengths of couplings, and a, 
/3, and 7 are the angles between the spins. 





FIG. 1: (Color online) (a) An isosceles triangle with spins 
(unit 3 vectors) at its vertexes. L and K are the strengths of 
couplings, and a, /3, and 7 are the angles between the spins. In 
the K < 0 case, there is no frustration and all the three spins 
are collinear. I11 the K > 0 case, there is a frustration. At 
fixed value of a = 2<f>, the minimum of the energy is attained 
when two other angles are equal to (f> if (b) L < 0 or to 7r — cj) 
if (c) L > 0. The spins are then coplanar. 


The energy of the system reads 

E = K cos a + L (cos (3 + cos7). (1) 

What are the values of a , /3, and 7 that minimize the en¬ 
ergy? The K < 0 case is trivial: the spins are collinear, 
because there is no frustration, that is the interaction 
energy can be minimized by pairs. If K > 0 (the case 
of frustration), then, at a fixed angle a = 2 <j>, the mini¬ 
mum value of the energy is attained when cos (3 + cos 7 
is maximum or minimum, depending on the sign of L. 
Since a + /3 + 7 < 27t, it is easy to see that, in the 
L < 0 case, the sum is maximum, if the spins are copla¬ 
nar and (3 = 7 = a/2 = <f> [Fig. 1(b)], and in the L > 0 
case, this sum is minimum, if the spins are coplanar and 
,0 = 7 = 7t — a/2 = 7 r — <j> [Fig. 1(c)]. 

The ground-state energy of the triangle equals to 

E gs = K cos2(t> — 2\L\cos(j). (2) 

The energy is minimum if cos <p = Thus, at |L| ^ 
2K, the ground state is determined by the angle (j) and 
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FIG. 2: (Color online) Ground-state configurations of a tri¬ 
angular ladder with (a) J < 0 and (b) J > 0. Jo > 0. The 
rotational angle along legs is equal to 2(j> = 2 arccos . All 
the spins are coplanar. Here, the spins are in the plane of 
the ladder but the polarization plane can be arbitrary, for in¬ 
stance, perpendicular to the ladder. The orientation of one 
arbitrarily chosen spin and the chirality (the same for both 
legs) can be arbitrary as well. 


the spins are coplanar (the polarization plane, the ori¬ 
entation of one of the three spins, and the chirality are 
arbitrary), otherwise the spins are collinear. 

Consider now a triangular ladder (Fig. 2). It is com¬ 
posed of identical triangle clusters. Each J-bond belongs 
to two such clusters; therefore, one should minimize the 
energy of a triangle cluster with the interactions Jo and 
J/2. In view of the previous result, we can assert that 
the ground state is spiral and determined by the angle 
c l ) = arccos jy-, if |J| < 4J 0 , otherwise all the spins are 
collinear. Two examples of the ground-state ordering for 
a triangular ladder are shown in Fig. 2. It should be 
emphasized that all the spins are coplanar but the polar¬ 
ization plane is arbitrary. Since the angle (j> depends on 
the strengths of couplings only, the spin ordering is, in 
general, incommensurate along legs. 

It is of interest to consider a 2D-lattice composed of 
two types of triangular ladders with nearest-neighbor in¬ 
teraction J 0 along legs and J\ and J 2 along rungs (Fig. 3). 
The lattice can be partitioned into identical two-triangle 
clusters. Each Ji- and J 2 -bond belongs to two such clus¬ 
ters; therefore, one should minimize the energy of a two- 
triangle cluster with the interactions Jo, Ji/2, and J 2 / 2 . 
The Jo < 0 case is trivial: there is no frustration. Con¬ 
sider the Jo > 0 case. Let the angle between the vec¬ 
tors at the ends of the Jo-bond be fixed and equal to 2<f 
(0 ^ </> ^ 7t/2). To minimize the energy, the angles be¬ 
tween neighboring spins joint by the Ji-bonds (z = 1, 2) 
should be equal to <f> if Ji < 0 or to 7r — </> if J t > 0. The 
energy per cluster is 

Egs = Jo COS 2(j) (| J11 -f- IJ2I) COS (f). (3) 



FIG. 3: (Color online) Left: A generalized triangular lattice 
and a spin configuration on it for Ji < 0 and Ji > 0 (Jo > 0, 
| Ji | + [ Ja | ^ 4Jo). A two-triangle cluster is shown by bold 
lines. Right: Square of spiral phases for the generalized tri¬ 
angular lattice. As the border of the square is approached, 
4> = arccos - tends to zero and the spiral phases con¬ 
tinuously evolve to their collinear limits. 


The angle <j> that minimizes this energy is given by 

I Ji I + \h\ 


COS0 = 


4J 0 


( 4 ) 


If | Ji | -F | J 2 1 < 4 Jo, then the spin ordering is spiral and 
determined by the angle 4> and by the signs of Ji and J 2 , 
otherwise all the spins are collinear. In accordance with 
the signs of Ji and J 2 , there are four spiral phases. The 
square of spiral phases is shown in Fig. 3 (right panel). 
One of these are depicted in Fig. 3 (left panel). As the 
border of the square is approached, <f> tends to zero and 
the spiral phases continuously evolve to their collinear 
limits where the spins belonging to the same leg have the 
same direction, although J 0 > 0. If J 2 = J\ = Jo > 0, 
then we have a well-known 120 °-phase on the triangular 
lattice. 

Let us pass to a 3D zigzag ladder lattice, shown in 
Fig. 4. We refer to it as a generalized hollandite lat¬ 
tice. The ground states of a classical Heisenberg model 
on the hollandite lattice (dark cyan and green rungs 
are identical) were studied in Ref. 21]. There one can 
find an overview of magnetic compounds with this lat¬ 
tice (see also Ref. [Hj]). Since the hollandite lattice is 
not a Bravais one, a straightforward application of the 
Luttinger-Tisza method is impossible; therefore, the au¬ 
thors used a modification of it to find the ground-state 
energy; they also performed numerical simulations. The 
cluster method makes it possible to solve this ground- 
state problem in a very simple and clear way. 

The generalized hollandite lattice is composed of three 
types of zigzag ladders. Let the strengths of coupling 
between neighboring spins along legs of all the ladders be 
Jo (blue bonds) and J\ (red bonds), J 2 (green bonds), 
and J 3 (dark cyan bonds) along their rungs (Fig. 4). It 
is easy to see that this lattice can be partitioned into 
identical three-triangle clusters shown in Fig. 4. 

Each Ji-, J 2 -, and J 3 -bond belongs to two such clus¬ 
ters; therefore, one should minimize the energy of a three- 
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FIG. 4: (Color online) The generalized hollandite lattice is 
composed of identical three-triangle clusters. 


triangle cluster with the interactions Jo, Ji/ 2 , J2/2, and 
J3/2. In the Jo < 0 case, there is no frustration. Con¬ 
sider the J 0 > 0 case. Let the angle between the spins 
at the ends of the Jo-bond be fixed and equal to 2 <j> 
(O^0^7r/2). To minimize the energy, the angles be¬ 
tween neighboring spins joint by the J^-bonds [i = 1 - 3 ) 
should be equal to 0 if Ji < 0 or to 7r — 0 if Ji > 0 (see 
Fig. 1 ). The energy per cluster (or per site) equals to 

E gs = J 0 cos 2 ^- (| Ji | + |J 2 | + |J 3 |)cos 0 . ( 5 ) 


The angle 0 that minimizes the energy is given by 


cos 0 


[Ji| + IJ2I + l^g] 

4 J 0 


( 6 ) 


It is clear that all the spins on the lattice should be 
coplanar that is parallel to a plane but this plane can be 
arbitrary. If the plane, the orientation of one arbitrarily 
chosen spin, and the chirality are fixed, then the distri¬ 
bution of angles completely determines the ground-state 
spin configuration of the lattice at fixed couplings. 

Each set of signs of the interactions Ji, J 2 , and J3 
(Jo > 0) corresponds to a spiral phase and its collinear 
limit. Since there are eight sets of signs, there are also 
eight spiral phases and eight collinear limits. The oc¬ 
tahedron of spiral phases and two spiral orderings with 
their collinear limits are shown in Fig. 5 . 

The cluster shown in Fig. 4 gives the ground-state or¬ 
derings of classical Heisenberg spins on many other zigzag 
ladder lattices (see Fig. 6). Eq. (6) and the octahedron 
of spiral phases (Fig. 5 ) are common for all these lat¬ 
tices. Hexagonal zigzag ladder lattices similar to those 
shown in Figs. 6(a) and 6(b]_have been identified in many 
magnetic compounds [ 13 h 2 C|. 

It is also interesting and instructive to consider a gener¬ 
alized triangular lattice composed of three types of zigzag 
ladder (Fig. 7 , left panel). This lattice cannot be parti¬ 
tioned into identical two-triangle clusters. Three types 
of two-triangle clusters have to be considered. It can 



FIG. 5: (Color online) Left: Octahedron of spiral phases for 
the classical Heisenberg model on the generalized hollandite 
lattice. Right: Two examples of spiral orderings and their 
collinear limits. Here, all the spins are parallel to the plane 
of the figure but this plane can be arbitrary. Cyan and black 
circles denote sites over and under the plane of the figure, re¬ 
spectively. Light gray arrows represent spins next to “black” 
ones along legs. 


be shown that a single-angle spiral ordering of classical 
Heisenberg spins is also possible in some regions of the 
parameter space of this model. Let us distribute the 
energy of each rung between two triangular plaquettes 
which share it in the manner shown in Fig. 7 (left panel). 
Coefficients cq for this distribution can be arbitrary but 
should satisfy the conditions 0 ^ cq ^ 1. 

The energies of two-triangle clusters can be written as 

El 2 = Jo COS 2012 — 2 («i| J1 1 + (1 — CK 2 )| J 2 |) cos 012, 

E23 = Jo COS 2023 - 2(a 2 |</2| + (1 - a 3 )|J 3 |) cos 0 23 , 

E31 = Jo cos 2031 — 2(a 3 | J 3 | + (1 — ai)| Ji|) cos 0 3 i. 

( 7 ) 

The condition for the minimum of these energies along 
with the condition 0i 2 = 0 23 = 0 3 i = 0 lead to the set 
of equations 


o?i | Ji | + (1 — a 2 )|J2| = ci^l^l + (1 — 0:3) | J 3 1 

„ iri , n „ mti l^il + 1*^21 + |> 73 1 roi 

= a 3 \J 3 \ + (1 - oq)| Ji| = -jj-• (8) 


The angle 0 that minimizes the energies is then given by 


COS0 


I Ji I + \'h\ + | J3I 

6J 0 


( 9 ) 


It it easy to show that the region in 
(Ji/- 7 o, J2/J0, J 3 /Jo)-space, where Eqs. (8) and the 
inequalities 0 ^ a* ^ 1 are satisfied, is the polyhedral 
cone determined by six vectors: (0, 2, 1), (0, 1, 2), (1, 0, 
2), (2, 0, 1), (2, 1, 0), and (1, 2, 0) in the first octant 
(Fig. 7 , right panel) and similar polyhedral cones in 
other octants (not shown in the figure). The single-angle 
spiral phases exist in the parts of this cones which 
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FIG. 6: (Color online) (a) - (d) Some hexagonal zigzag ladder 
lattices, (e) A zigzag ladder lattice based on the truncated tri- 
hexagonal tiling 4.6.12. The lattice (c) is a disordered mixture 
of the lattices (a) and (b). All these lattices are composed of 
identical clusters shown in Fig. 4. Cyan and black circles de¬ 
note sites over and under the plane of the figure, respectively. 
The legs of zigzag ladders are perpendicular to the plane of 
the figure. 




FIG. 7: (Color online) Left: A generalized triangular lattice, 
composed of three types of zigzag ladders, and three types of 
two-triangle clusters into which it can be partitioned. Right: 
Hexagonal pyramid (in red) of a single-angle spiral phase on 
this lattice for Ji is 0 {i = 0-3). 



FIG. 8: (Color online) A hexagonal zigzag ladder lattice real¬ 
ized in the S = § compound /3-CaCr2 04 and its partition into 
two types of three-triangle clusters. Cyan and black circles de¬ 
note sites over and under the plane of the figure, respectively. 
The legs of zigzag ladders are perpendicular to the plane of 
the figure. 


are bounded by the planes |Ji| + |J 21 + |J 3 I = 6 Jo- 
These are hexagonal pyramids, one of which is shown 
in Fig. 7 (right panel). Other parts of the polyhedral 
cones correspond to the collinear limits of the spiral 
phases. It seems that, out of the polyhedral cones, 
more complicated structures exist but here we study 
single-angle orderings only. It is easy to extend these 
results to arbitrary number of types of ladders. 

Finally, consider a hexagonal zigzag ladder lattice 
shown in Fig. 8 . This lattice is realized in the S = | 
compound /3-CaCr 2 04 18, 0 (however, the spin or¬ 
dering observed in this compound is more complicated 
than obtained here). In contrast to the lattices shown in 
Figs. 6 (a)- 6 (d), this lattice can be partitioned only into 
two types of three-triangle clusters (Fig. 8 , right panel). 
The ground-state energies for these clusters are given by 


Ei = J 0 cos20i - (2a| J 2 | + 2/31 J 3 | + \ J 4 1) cos(/>i, 

E 2 = Jo COS 2 (/) 2 

- [| Ji| + 2(1 — a)| J 2 I + 2(1 — /3)| J 3 |] cos02- (10) 


Minimizing these energies with respect to 0 1 = 02 = 0, 
we obtain 


COS0 = 


|Ji| + 2|J2| + 2 IJ 3 I + | J 4 1 
8 J 0 


( 11 ) 


The conditions 0 ^ a ^ 1 and 0 ^ /3 ^ 1 lead to the 
inequality 


||J4|-|Ji||<2(|J2| + |J 3 |). ( 12 ) 

This inequality along with the inequality | Ji | + 21J 2 1 + 
2|J 3 | + |J 4 | < 8 Jq determine the regions for existence of 
single-angle spiral phases. 

To summarize, we have rigorously proven that, in sys¬ 
tems of classical Heisenberg spins on zigzag (triangular) 
ladder lattices, there exist incommensurate ground-state 
spiral orderings (phases) characterized by a single angle 
and by the signs of the interactions between neighboring 
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spins (along ladder rungs). We propose to name these or¬ 
derings (phases) “single-angle spiral orderings (phases).” 
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